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We examine the quantum phase diagram of the fractional quantum Hall effect (FQHE) in the 
lowest two Landau levels in half-filled bilayer structures as a function of tunneling strength and 
layer separation, i.e., we revisit the lowest Landau level filling factor 1/2 bilayer problem and make 
predictions involving bilayers in the half- filled second Landau level (i.e., filling factor 5/2). Using nu- 
merical exact diagonalization we investigate the important question of whether this system supports 
a FQHE described by the non- Abelian Moore-Read Pfafhan state in the strong tunneling regime. 
In the lowest Landau level, we find that although in principle, increasing (decreasing) tunneling 
strength (layer separation) could lead to a transition from the Abelian two-component Halperin 331 
to non-Abelian one-component Moore-Read Pfafhan state, the FQHE excitation gap is relatively 
small in the lowest Landau level Pfafhan regime, and we establish that all so far observed FQHE 
states in half-filled lowest Landau level bilayers are most likely described by the Abelian Halperin 
331 state. In the second Landau level we make the prediction that bilayer structures would manifest 
two distinct branches of incompressible FQHE corresponding to the Abelian 331 state (at moder- 
ate to low tunneling and large layer separation) and the non-Abelian Moore-Read Pfaffian state 
(at large tunneling and small layer separation). The observation of these two FQHE branches and 
the possible quantum phase transition between them will be compelling evidence supporting the 
existence of the non-Abelian Moore- Read Pfaffian state in the second Landau level. We discuss our 
results in the context of existing experiments and theoretical works. 

PACS numbers: 73.43.-f, 71.10.Pm 



I. INTRODUCTION 



Two important developments have rekindled interest 
in the phenomena of even-denominator incompressible 
fractional quantum Hall states in two-dimensional high- 
mobility bilayer semiconductor structures. The first is 
the recent intriguing experimental observation by Luh- 
man et al} of two distinct even-denominator fractional 
quantum Hall effect (FQHE) states at filling factors 

= 1/2 and 1/4 in a very wide (~ 600 A) single quan- 
tum well (wide-quantum-well) structure at very high 
(> 40T) magnetic fields. The second development, moti- 
vated by implications for fault-tolerant topological quan- 
tum computation^'^, is the great deal of recent theo- 
retical and experimental interest in the possible non- 
Abelian nature of the v = b/2 second Landau level (SLL) 
FQHE state first observed by Willett et al.^ in the single- 
layer system in 1987 with subsequent confirming obser- 
vation8^~— over the years since. The FQHE at filling 
factor V = b/2 presents an amazing confiuence of ideas 
from condensed matter physics, conformal field theory, 
topology, and quantum computation.^ In particular, the 
fundamental nature of the experimentally observed 5/2 
FQHE, whether an exotic spin-polarized non-Abelian in- 
compressible paired state or a more common Abelian in- 
compressible paired state, has remained an enigma for 
more than 20 years. Although theoretical and numeri- 
cal work indicates that the 5/2 FQHE belongs to a non- 
Abelian Moore- Read Pfaffian^^ universality class, there is 



scant experimental evidence supporting this conclusioni^ 
(see Ref. 3 for a comprehensive review of non-Abelian 
physics and topological quantum computation). 

These two developments lead to important and inter- 
esting questions. One question is whether a non-Abelian 
V = 1/2 FQHE, i.e., the analog of the possibly non- 
Abelian = 5/2 = 2 + 1/2 SLL single-layer state^^, can 
exist in the lowest Landau level (LLL) under experimen- 
tally observable conditions. This question has a long his- 
|-Q^yi5 18 ^YiQ theoretical literature going back to the 
early 1990s and recent FQHE experiments at v = 1/2 
make it imperative that a theoretical analysis be carried 
out to achieve a proper qualitative understanding of cur- 
rent experiments 

Another important question concerns the non-Abelian- 
ness (or not) of the experimental 5/2 state which is of 
profound importance beyond quantum Hall physics^^^. 
Hence, it is useful to contemplate novel situations where 
the nature of the 5/2 state will manifest itself in a dra- 
matic, but hitherto unexplored, manner. The current 
work addresses the dichotomy (i.e., the 1/2 FQHE be- 
ing an Abelian Halperin 331 bilayer state^S versus the 
5/2 FQHE being a non-Abehan Moore-Read Pfaffiani^^^ 
single-layer state) between single-layer SLL physics ver- 
sus bilayer LLL physics by studying the bilayer FQHE 
in both the lowest and second LLs as a function of layer 
separation and tunneling. We ask whether bilayers can 
support both Pf and 331 FQHE. 

The proposed non-Abelian Moore- Read Pfaffian (Pf) 
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state is a weak-pairing single-layer (or one-component) 
FQHE state for half-filled LLs which, in principle, applies 
to any orbital LL (i.e., LLL as well as SLL). Thus, as a 
matter of principle a = 1/2 LLL single-layer Pf FQHE 
is certainly a possibility^^"— '^^^^ although it has never 
been observed experimentally. Theoretically, due to the 
differences in the electron-electron Coulomb interaction 
pseudo-potentials in the LLL compared to the SLL the 
two situations (i.e., 1/2 and 5/2) are quantitatively very 
different, and it is possible for the Pf FQHE to exist in 
the SLL, but not in the LLL (and vice versa). The best 
existing numerical woi]^^^^^ indicates that either the 
single-layer u = 1/2 LLL Pf state does not exist in na- 
ture or if it exists, does so only in rather thick 2D layers 
with an extremely small FQHE excitation gap, making 
it impossible or very difficult to observe experimentally. 
By contrast, the single-layer u = 5/2 SLL FQHE is ob- 
served routinely, albeit at low temperatures (^ lOOmiC), 
in high mobility (^ 10^ cm^ /Vs) samples, and with a 
rather small (but experimentally accessible) activation 
gap (~ 100-500mi^). In fact, it has been pointed out 
that the experimental u = 5/2 FQHE is always among 
the strongest observed FQHE states in the SLL. 

Instead of studying a single-layer 2D system, we con- 
centrate on the spin-polarized and density balanced bi- 
layer system assuming an arbitrary tunneling strength t 
and an arbitrary layer separation d in the lowest and 
second LLs. We confine our study to spin-polarized 
systems because all experimentally realized half-filled 
FQHE states appear to be fully spin-polarized, consis- 
tent with the theoretical expectation^^^^^-whether bi- 
layer two-component FQHE states at z^ = 1/2 or one- 
component single-layer FQHE states at u = 5/2. By 
density "balanced" we mean that each layer in the bilayer 
system has the same number of electrons. Furthermore, 
it should be noted that the tunneling strength t is propor- 
tional to the symmetric-antisymmetric energy gap. The 
experimentally realizable system we have in mind could 
either be a true bilayer double-quantum-well structure 
or a single wide-quantum-well (WQW) which manifests 
effective bilayer behavior where the self-consistent field 
from the electrons produces effective two-component be- 
haviori^ (see Fig. [H] in Sec. lVlIIl) . 

We numerically obtain the approximate quantum 
phase diagram for a 1/2 filled bilayer system in either 
the lowest or second Landau level in the t-d space us- 
ing the spherical system finite size exact diagonalization 
(Lanczos) technique, concentrating entirely on the Pf and 
the 331 FQHE phases. 

For the LLL, we revisit the u = 1/2 bilayer FQHE 
and carry out an extensive comparison with all existing 
bilayer = 1/2 FQHE experimental observations to as- 
certain any hint of the existence of a non-Abelian Moore- 
Read Pfaffian state for large values of t. No strictly 
single-layer system, e.g., a heterostructure or a not-too- 
thick quantum well, has ever demonstrated an incom- 
pressible FQHE at u = 1/2, instead manifesting only the 
compressible composite fermion2£i22 Fermi sea^^^^^. It is 



intuitive that our model system is an effective bilayer, or 
single-layer, u = 1/2 system for small, or large, values of 
t, and therefore by studying the quantum phase diagram 
as a function of t and d we hope to shed light on the 
possible existence of a single-layer = 1/2 FQHE in real 
systems. 

One of the results of the current work is that for the 
lowest Landau level bilayer FQHE system, (i) the re- 
cently observed WQW z^ = 1/2 FQHE^ii^ are strong- 
pairing Abelian Halperin 331 FQHE states^^'^^ which, 
however, sit close to the boundary between the Abelian 
331 and the weak-pairing non-Abelian Pfaffian FQHE 
state^^, and (ii) it may be conceivable, as a matter of 
principle, to realize the LLL u = 1/2 Pfaffian non- 
Abelian FQHE in very thick bilayers, but as a matter of 
practice, this is unlikely since the z^ = 1/2 FQHE gap is 
extremely small (perhaps zero) in the parameter regime 
where the Pf is more stable than the 331 phase. Our find- 
ings about the fragility of the LLL z^ = 1/2 non-Abelian 
Pf state are consistent with recent conclusion a^^i^^i^^i^^ , 
but our main focus in the current work is in understand- 
ing the z^ = 1/2 bilayer quantum phase diagram treating 
tunneling layer separation and individual layer width 
w of the 2D system as independent tuning parameters of 
the Hamiltonian. 

For the second Landau level, we predict that bilayers 
can support both the non-Abelian Moore-Read Pfaffian 
and the Halperin Abelian 331 FQHE and that there could 
be a novel quantum phase transition, both as a function 
of tunneling strength (at constant layer separation) and 
of layer separation (at constant tunneling), between the 
two-component 331 Abelian and the one-component Pf 
non-Abelian states in a bilayer SLL system. We show 
that tuning the inter-layer tunneling and/or layer sep- 
aration would lead to a transition between the Abelian 
and the non-Abelian SLL FQHE, which should be observ- 
able experimentally in standard FQHE transport experi- 
ments. In particular, we predict that in realistic systems 
with finite single-layer width, the SLL bilayer state would 
manifest two distinct FQHE phases separated by a region 
of finite inter-layer separation and tunneling. Existence 
of two distinct incompressible FQHE bilayer states at to- 
tal filling factor v = 5/2, connected possibly by a quan- 
tum phase transition, is a clear (experimentally testable) 
prediction of our theory. The observation of such a quan- 
tum phase transition, originally predicted as possible by 
Read and Green^^ (Fig. 1 in Ref. 34) under general the- 
oretical considerations (but never before demonstrated 
to be feasible under realistic conditions), would strongly 
suggest the existence of a non-Abelian 5/2 state since 
the two distinct FQHE phases in the same sample both 
cannot conceivably be Abelian 331 states. 

We first present a background for our work in Sec- 
tion ini and then describe our theoretical model in terms 
of a Hamiltonian and introduce all relevant parameters 
in Section IIIII Next we revisit the lowest Landau level 
problem in Section HVl before tackling the second Landau 
level problem in Section|Vl since the LLL problem, due to 
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its long history, is easier to understand and will provide a 
proper context and atmosphere when discussing our SLL 
results. In Section ITV CI we discuss in detail how we con- 
nect our results with current and previous experimental 
FQHE bilayer LLL results. Furthermore, in Section IV Dl 
we discuss am important issue regarding bilayer FQHE 
systems in higher Landau levels-a difficulty or ambiguity 
that is quite subtle and has not been discussed previously 
as far as we know. Lastly we present our conclusions in 
Section ED 



II. BACKGROUND 

Before describing and presenting our work, it is useful 
to provide a brief background of bilayer FQHE, a subject 
with a long history, in order to set a context for our work. 

The two candidate wavefunctions we consider and com- 
pare in this work, with respect to = 1/2 and 5/2 bilayer 
incompressible states, the Halperin 331 state^^^^^ and 
the Moore-Read Pfaffian statei^ were proposed in 1983 
and 1991, respectively. The 331 wavefunction is a two- 
component strongly paired Abelian state of two Laugh- 
lin^^ phases in each layer whereas the Pfaffian is a one- 
component weak pairing non-Abelian superconducting 
state of chiral p-wave symmetry. Both of these states are 
allowed incompressible fractional quantum Hall (FQH) 
phases at half-filled Landau levels which for our purpose 
could be either v = 1/2 or 5/2(= 2 + 1/2). The quali- 
tative difference between these two incompressible FQH 
states is that the 331 is an Abelian two-component state 
whereas the Pfaffian is a non-Abelian one-component 
state. 

It was first explicitly shown by Yoshioka, Girvin, and 
MacDonald^^'^^ that a bilayer system at (total) half- 
filHng could support a 331 incompressible state for a layer 
separation d ~ / (where / is a characteristic length scale 
called the magnetic length, defined below), i.e., when 
intra- and inter-layer correlations are comparable. Later, 
He et al}^'^^ carried out a detailed quantitative analysis 
of the possible existence of 331 FQHE in z^ = 1/2 bilayer 
structures, making specific predictions for layer separa- 
tion (d) and layer width (w) values where experimentally 
observable incompressible states may exist. This lead to 
the observation of u = 1/2 bilayer 331 FQHE by Eisen- 
stein et al.^^. Parallel to the Eisenstein et at. obser- 
vation of the 331 FQHE in bilayers, Shayegan et a/.^^^^^ 
observed a well-defined v = 1/2 FQHE in wide-quantum- 
wells. This WQW observation oiu = 1/2 FQHE was at 
first attributed to the one-component Moore-Read Pfaf- 
fian state by Greiter, Wen, and Wilczek^^, who argued 
that the occurance of a z^ = 1/2 FQHE in a single quan- 
tum well, rather than in double quantum wells, implies 
a one-component rather than two-component nature of 
the underlying incompressible state. Later theoretical 
work by He et al}^ and further experimental work by 
Shayegan et al.^^'"^^ decisively established that the wide- 
well V = 1/2 FQHE is, in fact, a manifestation of the 



two-component 331 rather than the one-component Pfaf- 
fian state. The existence of a two-component 331 FQHE 
in wide single wells becomes possible by virtue of the 
self-consistent Hartree electric field arising from the elec- 
trons themselves which, for sufficiently large quantum 
well width {W) and carrier density (n), could lead to the 
single wide-quantum-well acting effectively as two dis- 
tinct electron layers localized near the well boundaries 
(with density n/2 each) with a potential barrier sepa- 
rating them in the middle. Such a two-component 331 
FQHE description of observed v = 1/2 incompressibility 
has become well-accepted, and the fact that no v = 1/2 
FQHE has ever been observed in relatively thin single 
quantum well structures or in single heterostructures has 
further reinforced the idea that the = 1/2 FQHE is 
due to the formation of the two-component Halperin 331 
state in wide-quantum-wells^. 

From the perspective oi v = 1/2 FQHE, wide- 
quantum-wells should be considered as two-component 
bilayer systems with non-zero interlayer tunneling. For 
relatively weak, or strong, tunneling between the lay- 
ers, the system behaves as a two, or one, component 
system, and the question of the existence of Abelian 
331 or non-Abelian Pfaffian v = 1/2 FQHE then, in 
some sense, boils down to the existence of incompress- 
ible FQHE in the weak or strong tunneling limit. Un- 
fortunately, how strong an interlayer tunneling is strong 
enough to render the system into a one-component Pfaf- 
fian FQHE is a quantitative question, which can only 
be addressed through detailed numerical study. Such a 
numerical study is the main goal of this paper. In the 
process, we also ask the question of one-component ver- 
sus two-component (i.e., Pfaffian versus 331) FQHE in 
bilayer v = h/2 case, where the existence of the one- 
component Pfaffian v = h/2 FQHE is reasonably well- 
established theoretically^^'^''^'^^^^, although actively de- 
bated^^^^. The bilayer v = b/2 case was never stud- 
ied before in the literature whereas there were only two 
studies of the bilayer z^ = 1/2 case comparing 331 versus 
Pfaffian state. The first is by He et al}^ and the second 
is by Nomura and Yoshioka^^. Our work for v = 5/2 
is distinct, and our work for v = 1/2 transcends that 
of the earlier work in being much more complete. Very 
recently, Papic et al?^ investigated this problem using a 
somewhat different model. 

As mentioned in the Introduction (Sec. [Ij), our work 
is partially motivated by the recent experimental obser- 
vation of Luhman et al} who found a z^ = 1/2 FQHE 
in side single wells with relatively strong tunneling (i.e., 
large symmetric- antisymmetric energy splitting). An in- 
teresting question is whether the v = 1/2 FQHE observed 
in the Luhman et al. experiment is a two-component 331 
state or a one-component Pfaffian state. We do not study 
the z^ = 1/4 FQHE observed by Luhman et al. which has 
recently been discussed by Papic et al.'^^ 

We mention that all our theoretical work assumes com- 
plete spin-polarization of the electrons and considers only 
the balanced case where the average electron density is 
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the same in each layer. We also neglect all Landau level 
coupling effects, and as such our work does not distin- 
guish between the non-Abelian Pfafhan and the non- 
Abehan anti-Pfaffian^^^^ states at = 1/2 or 5/2. (Note 
that the neglect of LL mixing effects may not be a par- 
ticularly good assumption for the 5/2 FQHEi^i^.) The 
reason for our considering full spin-polarization is that 
prior theoretical work by Morf^^ and by Feiguin et al?^ 
indicates that the i/ = l/2or5/2 state is likely to be fully 
spin-polarized. If the bilayer incompressible states turn 
our to be unpolarized or partially polarized, our work 
simply would not be valid. 

With this background, we study the stability of the 
331 and the Pfaffian state in v = 1/2 and 5/2 bilayers 
in the presence of finite interlayer tunneling (t), inter- 
layer separation (d), and layer width (w). Our goal is 
to obtain an appropriate zeroth-order quantum phase di- 
agram in the t-d-w space for u = 1/2 and 5/2 bilayer 



FQHE. Given that we calculate and compare the exact 
many-body ground state for many individual values of t, 
and we are forced to use a rather modest system 
size for A/" = 8 electrons (i.e., 4 in each layer) in all our 
exact diagonalization work. Past experience shows that 
a 8-electron system is quit adequate for qualitative un- 
derstanding as long as precise thermodynamic values of 
excitation gap or ground state energy are not desired. 



III. THEORETICAL MODEL 

We use the simplest model Hamiltonian H incorporat- 
ing both finite tunneling t and finite layer separation d 
(as well as finite layer width w-oi finite thickness) for 
our bilayer FQHE system: 



N 
i<j 

Ti and Ti are the position of the i-th electron in the right and left layer, respectively. In Eq. ([T]), l^ntra(^) = 
e^/ {K\/r'^ + w'^) (we use the Zhang-Das Sarma^*" potential to model the single layer quasi-2D interaction) and 
Mnter(^) = 6^/ (/^V^^ + d^) are the intralayer and interlayer Coulomb interaction incorporating a finite layer width 
w and a center-to-center interlayer separation d {> w hy definition). The x-component of the pseudo-spin operator 
{Sx)\a.yev (written in the layer basis representation) controls the tunneling between the two quantum wells with large 
t denoting strong tunnehng. 

Note that in the bilayer problem when the electron density is balanced in each layer, i.e., total number of particles 
in each layer is A^/2, there are essentially two natural Hilbert space representations. The layer basis (in which Eq. [T]is 
written) or the symmetric- antisymmetric basis where CmS = {cmR + CmL) / a/2 and CmA = {cmR — CmL) / a/2 destroy an 
electron in the symmetric {S) and antisymmetric {A) superposition states, respectively, where m is angular momentum 
and CmR and CmL destroy an electron in the right and left quantum well, respectively. S{A) can be considered to 
be an effective pseudo-spin index for the bilayer system. In the symmetric- antisymmetric basis, the Hamiltonian is 
written as 

^ =\ (miai,m2a2|V'|m4Cr4,m3Cr3)c]^^^^C^^,C^4a4Crn3a3 - \ Y^{cLs^mS - C^mA^^mA) • (2) 

{mi,cri=A,S} rn 



As written above in Eq. [H the intra-layer Coulomb po- 
tential energy between two electrons is F = l^ntra and 
the inter-layer potential energy is V = l^nter- A differ- 
ence between the two representations is that the tunnel- 
ing operator in the layer-basis can be written as (-S^c) layer 
while, in the symmetric- antisymmetric basis, the tunnel- 
ing operator is {Sz)sAS \ Y^m^^Ls^rnS - CmA^mA) 
^Em(cLi?CmL ^C^mL^^mR) = (^x)layer. The layer- and 

symmetric-antisymmetric-bases are related through a 
pseudo-spin rotation. Of course, the choice is a matter 
of personal taste and convenience and one should really 
appeal to a physical explanation to understand the tun- 
neling: no matter the basis choice, the tunneling term 



controls the probability that electrons jump back and 
forth between the two layers keeping the total number 
of electrons in each layer fixed at A'/2, i.e., keeping the 
bilayer system density balanced. 

We numerically diagonalize H for finite N assuming 
specific values of and t (each expressed through- 

out in dimensionless units using the magnetic length 
/ = {ch/eBy^'^ as the length unit and the Coulomb en- 
ergy e^/{nl)^ where n is the background dielectric con- 
stant of the host semiconductor, as the energy unit; e is 
electron charge, c is the speed of light in vacuum, and B 
is the magnetic field strength). We utilize the spherical 
geometry where the electrons are confined to a spherical 
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surface of radius R = N^j) is the total magnetic 

flux piercing the surface (A^^ is an integer according to 
Dirac), and the filling factor in the partially occupied 
Landau level (whether LLL or SLL) is liniAr^oo • 
To fully consider tunneling, the Hilbert space must con- 
tain basis states with TV and N — N electrons in the right 
and left layers for TV G [0, TV]. 

Following the standard well-tested proce- 
dure^"— i^i^i^*^ used extensively in the FQHE 
literature, we calculate the overlap between the exact 
numerical TV electron ground state wavefunction of the 
Coulomb Hamiltonian defined by Eq. ^ and the candi- 
date TV electron variational states which are the Abelian 
Halperin 331 strong-pairing^^ and the non- Abelian 
Moore-Read Pfaffian weak-pairing^^ wavefunctions: 



N/2 N/2 N 

M/331 = ll{'i-zjfl[Czi-zjfl[iz,-zj) (3) 

i<j i<j i,j 




respectively, where z = x — iy is the electron coordinate 
in the x-y plane. Intuitively, the "331" in ^331 origi- 
nates from the exponents for {zi — Zj)^ {zi — ij), and 
{zi — Zj)^ respectively, i.e., 3, 3, and 1. (Note that other 
Halperin 331-type wavefunctions^^^^^ have been consid- 
ered for other filling factors, e.g., total v = 1^ and re- 
cently the total i/ = 5/2 + 5/2 = 5 bilayer FQHE has 
been considered^^, however, none of these states directly 
apply to the situation that we are investigating of a 
half-filled LLL or SLL.) Both ^331 and ^pf are written 
above in the layer representation, where the 331 state 
pairs electrons between layers and pairs electrons 
in a single-layer, since writing the wavefunctions in real 
space and appeahng to the layer basis representation is 
much more intuitive. However, we are working within the 
"balanced" density situation so the Pf wavefunction, in 
particular, should be thought of as a wavefunction that 
describes pairing among electrons in the symmetric state 
{{cmR^ CmL) / V^) siucc the clcctrous are never physically 
completely in a single layer, right or left. 

We emphasize that we only consider the above two can- 
didate wavefunctions. It is well known^^'^^ that a com- 
posite fermion liquid phase (or composite fermion Fermi 
sea) exists ior v = 1/2 one-component systems both from 
extensive theoretical analysis and experimental observa- 
tions^. On the other hand, it is strongly suspected that 
for one-component systems at = 5/2, the ground state 
is the Moore-Read Pfaffian phase. Again, this is known 
from theoretical and experimental works— i^ii^ii^iiiii^. 
There are other candidate wavefunctions that can be 
written down besides ^331 and ^pf. Namely, one could 
consider ^222 and ^440 (both obvious generalizations of 
naming convention used for ^331) that describe a pseudo- 
spin unpolarized composite fermion Fermi sea and two 
uncorrelated one-quarter-filled composite fermion Fermi 
seas, respectively^^. However, those states have different 



"shifts" (see below) in the spherical geometry than one 
another and ^331 and ^pf and cannot be compared on 
an equal footing^^ and comparing them each to the exact 
state ^0 would require many additional calculations that 
are simply beyond the scope of this work. 

After diagonalizing H we ensure that (i) the ground 
state is homogeneous, i.e., has total orbital angular mo- 
mentum L = and is therefore an incompressible state, 
and (ii) there is a gap, the FQHE excitation gap, sep- 
arating the ground state from all excited states. We 
are concerned with the Pf and 331 variational states 
for total filling factor 1/2 in either the lowest or sec- 
ond Landau level, so TV^ = 2TV — 3, the "-3" is known 
as the "shift" and is a consequence of the curvature of 
the spherical geometry in which we work. Throughout, 
we consider (TV, TV^) = (8,13) and mention that due to 
the pseudo-spin component present in the bilayer prob- 
lem the Hilbert space is large (more than 10^ states for 
TV = 8 and over 7 x 10^ for TV = 10). Usually when 
one is utilizing exact diagonalization one wishes to con- 
sider many different system sizes and extrapolate the fi- 
nite size results to the infinite system size through some 
sort of finite-size scaling. The computational difficulty 
of our problem does not allow this. TV = 6 electrons at 
TV^ = 2TV — 3 is aliased with a one-component FQHE 
corresponding to = 2/3-remember the filling factor 
V = limAT^oo N/Nfjy^ so for certain combinations of TV and 
TV^ there can be correspondence between the finite sys- 
tems and two distinct filling factors (the so-called aliasing 
problem). Hence, the TV = 6 case will produce ambigu- 
ous results (and, in fact, this aspect raises some ques- 
tions about the work of Nomura and Yoshioka^^ who 
used TV = 6 to do a similar investigation for v = 1/2 
bilayers). TV = 10 is too big for us to exactly diagonalize 
(mentioned above) so we are left with only being able 
to consider TV = 8 and no finite size scaling is possible. 
Although this is a drawback, it is not uncommon in the- 
oretical FQHE studies to use a modest system size, but 
many different sets of physical parameters (i.e., t, (i, w for 
our case) to bring out qualitative features. Furthermore, 
since the theoretical techniques are standard, we do not 
give the details, concentrating instead on the results and 
their implications for bilayer FQHE experiments. 

We investigate this system with the usual probes used 
in theoretical FQHE studies by calculating: (i) wavefunc- 
tion overlap between variational ansatz (^pf and ^331 
states) and the exact ground state ^0 of iJ-an overlap 
of unity or zero indicates that the physics is or is not 
described by the ansatz; (ii) expectation value of the ex- 
act ground state of (TV^ — TV^)/2 , where TV^ and TV^i 
are the expectation values of the number of electrons in 
the symmetric or the antisymmetric states, respectively- 
a value of zero or TV/2 indicating the ground state to be 
two- or one-component; and (iii) energy gap (provided 
the ground state is a uniform state with L = 0)-a non- 
zero excitation gap indicating a possible FQHE state. In 
other words, we ask: (i) What is the physics (i.e., 331 
or Pf)?; (ii) Is the system one- or two-component?; (iii) 
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Will the system display FQHE (i.e., is the system com- 
pressible or incompressible)? 

The calculated overlap and gap determine the nature 
of the FQHE and its strength in our theory. We opera- 
tionally define the system to be in the 331 or Pf phase 
depending on whether the overlap between the exact 
ground state and the 331 or Pf wavefunctions is larger, 
i.e., if (^o|^33i) > (^ol^Pf) then the system is said to be 
in the 331 phase and vice versa. We emphasize that our 
work is a comparison between these two incompressible 
states only, and we cannot comment on the possibility 
of some other state (i.e., neither 331 nor Pf) being the 
ground state. We do, however, believe that if the system 
is incompressible at a particular set of parameter values 
(i.e., (i, t, etc.), it is very likely to be described by one 
of these two candidate states, 331 or Pf. We cannot, how- 
ever, rule out the possibility that the real system has a 
compressible ground state (without manifesting FQHE), 
e.g., a composite fermion Fermi sea, not considered in 
our calculation. This is more likely to happen when our 
calculated excitation gap is very small. 



IV. LOWEST LANDAU LEVEL RESULTS 
A. What is the physics?— lowest Landau level 

The calculated wavefunction overlap between the exact 
ground state and the two appropriate candidate varia- 
tional wavefunctions (^pf and ^331) as a function of dis- 
tance d and tunneling energy t is shown in Fig.[TJ First we 
focus on the situation with zero width w = (Fig. [Il^a)) 
and concentrate on the overlap with ^pf. In the limit of 
zero tunneling and zero separation, the overlap with ^pf 
is zero and quickly jumps to approximately 0.9 when the 
tunneling is increased to only a small amount of approx- 
imately 0.025, and for increasing tunneling, (^ol^Pf) re- 
mains relatively constant. In fact, this zero d and large 
tunneling result should be compared to the single layer 
results we have given earlier (Fig. 9(a) in Ref. 25). For 
(i ^ 0, in the weak tunneling limit the overlap with ^pf 
remains very small. However, in the large tunneling limit, 
increasing d only reduces the overlap marginally and the 
Pf description remains quite good even for large d. This 
shows that the strong-tunneling (one-component) regime 
is well described by the Pf state. It should be noted, how- 
ever, that the overlap between ^0 and ^pf is never much 
above 0.9-this should be compared to the SLL where 
we know2^5^i2iii^ii^5iii^ that the Pf wavefunction is a 
good physical description in the small layer separation 
and large tunneling limit (see Section IV Ap . 

Next we consider the overlap between ^0 and ^331. 
In the zero tunneling limit, as a function of we find 
the overlap starts very small, increases to a moderate 
maximum of ~ 0.80 at d ~ 1 before achieving an essen- 
tially constant value of ~ 0.6. For d > 4 the overlap in- 
creases with increasing tunneling to a maximum of nearly 
0.8. Thus, the weak-tunneling (two-component) regime 



is well described by ^331 and the 331 state is very robust 
to tunneling when layer separation is large. For small 
layer separation {d < 2), non-zero tunneling very quickly 
suppresses the overlap between ^331 and the exact state 
to a very small value. 

We know from earlier work^^i^ that finite layer width 
within a single quantum well, at least in the large tun- 
neling and small layer separation limit, i.e., the one- 
component limit, may enhance the overlap between ^0 
and ^pf. In Fig. [2](a) and (b) we show (^o|^33i) and 
(^ol^Pf), respectively, versus layer separation d and sin- 
gle layer width w for small tunneling t = and large tun- 
neling t = 0.2. For the small tunneling situation where 
331 is a good ansatz, we see clearly that for w = the 
overlap is maximum for a finite value of d 2. How- 
ever, for finite w we find that the position of maximum 
overlap does not change much and, in fact, for large w 
the maximum overlap obtains essentially for the w ^ d 
condition. 

In the case of ^pf (Fig. [2](b)) we find a result similar 
to the single- layer finite-thickness results (cf. Ref. |25[ ) 
where the overlap with the exact state is nearly constant 
for increasing w and d. For w = dwe expect this behavior 
but it is interesting to note how robust the overlap is for 
finite w and large d. 

In Fig. [Hb)-(d) we consider the overlap between the 
exact state and the 331 and Pf states as function of sep- 
aration and tunneling for finite layer widths of w = 0.6, 
1.2, and 2.4. If we consider first (^ol^Pf) we see that 
not much changes for different values of w as one would 
expect from examining Fig. [2fb). The only real qualita- 
tive change is that for w = d for w — 0.6 and 2.4 the 
overlap is non-zero, albeit very small, compared to the 
w — ^ case where it is identically zero. However, for 
iL' = 1.2 = d the overlap is zero, thus, we suppose that 
the small but non-zero overlap for d = is simply a fi- 
nite size effect. For the overlap between the exact state 
and the 331 variational state we find that finite single 
layer width has two effects. First, it slightly lowers the 
maximum overlap obtained (this is evident in Fig. [2]) and 
reduces the area in d-t phase space where the 331 state is 
a good ansatz. In essence, for finite w the Pfaffian phase 
pushes out the 331 phase, at least as far as wavefunction 
overlap determines the phase. 



B. Is the system one- or two-component?— lowest 
Landau level 

Due to the potential confusion between whether the 
pseudo-spin operator {Sz)sas or (iSa^) layer controls the 
tunneling we will report our results in terms of (Ns — 
Na)/2 where Ns and Na are the expectation values of 
the total number of particles in the symmetric S or an- 
tisymmetric A state, respectively. This is the more pre- 
cise variable since whether we are calculating {Sz)sas or 
{Sx)\asyev the answer is always equal to {Ns — Na)/2. 

Our conclusion so far, based on overlap calculations. 



7 



w =0, LLL w =0.6, LLL 




FIG. 1: (Color online) Wavefunction overlap between the exact ground state, ^o, and the Pf state, ^pf (blue (dark gray)), and 
331 state, ^331 (yellow (light gray)), as a function of layer separation d and tunneling amplitude t for the = 1/2 lowest LL 
with N = 8 electrons and single layer width (a) w = 0, (b) w = 0.6, (c) w = 1.2, and (d) w = 2.4 (d > w necessarily). 




FIG. 2: (Color online) Wavefunction overlap between the exact ground state ^0 and ^pf at strong tunneUng t = 0.2 (a) and 
^331 at zero tunneling t = (b) as a function of separation d and single- layer well width w (where d > w) for the z/ = 1/2 
lowest LL with N = 8 electrons. White corresponds to an overlap of unity while black corresponds to an overlap of zero. The 
dashed red line is the condition w = d and for w > d the bilayer system is undefined, i.e., the single- layer width cannot be 
larger than the layer separation. 



is completely consistent with the calculated expectation 
value of {Ns — Na)/'2 which is essentially the order pa- 
rameter describing the one-component to two-component 
transition, i.e., {Ns — Na)/2 ^ N/2 describes a one- 
component phase since in that case N = TV^, whereas 



{Ns — Na)/2 ^ describes a two-component phase since 
in that case Ns = Na = N/2. Note that in our TV = 8 
system the maximum and minimum value of the pseudo- 
spin expectation value is 4 and 0, respectively. Fig.[3l^a)- 
(d) shows the calculated {Ns — Na)/2 as a function of 
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w =0, LLL w =0.6, LLL 




(c) (d) 



FIG. 3: (Color online) [Ns — Na)/2 (or {{Sz)sas) when the Hamiltonian H is written in the SAS-basis) as a function of layer 
separation d and tunneling amplitude t for = 1/2 and N = 8 electrons for single-layer width (a) w = 0, (b) w = 0.6, (c) 
w = 1.2 and w = 2.4. 



d and t for w = 0, w = 0.6, w = 1.2, and w = 2A . 
First we consider w = (Fig. [3l^a)). For d = and t = 
the system is two-component and the system is SU(2) 
symmetric. Only a small amount of tunneling is required 
to very quickly push the system to be one-component. 
As d is increased, more tunneling is required to make the 
system one-component. Of course, this can be readily un- 
derstood physically: for non-zero layer separation d the 
electrons in the symmetric state pay a higher Coulomb 
energy price than electrons in the antisymmetric state 
since 1/r > l/^/r^T~d?^ for d ^ 0. 

For non-zero single-layer width (finite w) the general 
behavior outlined above does not change appreciably, cf. 
Figs.[3l^b)-(d). The only real difference is that for increas- 
ing less tunneling is required, at the same value of layer 
separation d, to make the system one-component. That 
is, it is easier to push the system into one-component 
behavior via tunneling when w ^ 0. 

Comparing these results to our overlap results shown 
in Fig. [T] we see that when the system is effectively two- 
component for weak tunneling the Halperin 331 state has 
a higher overlap with the exact state than the Moore- 
Read Pfaffian state. This is expected since the Pf state 
describes a one-component state. On the other hand, 
when the system is effectively one-component for strong 
tunneling the one-component Pf state has a higher over- 
lap with the exact state than ^331 does. Looking 
closer, one finds that when {Ns — Na)/'^ is increased 
from zero towards {Ns — Na)/2 ~ 2.5 the overlap with 



the exact state switches from being higher with the two- 
component 331 state to the one-component Pf state (this 
is true irrespective of w). In other words, the two- 
component description of the exact state provided by 
^331 survives until nearly 60% more electrons are in the 
symmetric state than the antisymmetric state-the two- 
component 331 state is surprisingly robust to tunneling 
that drives the system towards one-component behavior. 



C. Will the system display the FQHE?-lowest 
Landau level 

Wavefunction overlap and pseudo-spin are only two 
properties that elucidate the physics. Another property 
is the energy gap above the L = ground state in the ex- 
citation spectra; a crucial characteristic that determines 
the incompressibility or the robustness of the FQHE. In 
Fig. Ill^a)-(d) we show the energy gap, defined as the dif- 
ference between the first excited and ground state en- 
ergies at constant A^^, as a function of d and t for the 
v = 1/2 lowest LL system for (a) w = 0^ (h) w = 0.6, (c) 
w = 1.2, and (d) w = 2.4. It is clear that for finite d and 
t there is a 1/2 FQHE with a finite gap. At the SU(2) 
symmetric point the energy gap is vanishingly small. In- 
terestingly, the energy gap has a peak in t-d space-a ridge 
along which the energy gap is maximum. Generally, as 
w is increased the value of the energy gap decreases as is 
expected since finite width of a single quantum well re- 
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FIG. 4: (Color online) FQHE energy gap as a function of layer separation d and tunneling amplitude t for the = 1/2 lowest 
LL system with N = 8 electrons. We consider single quantum well widths of (a) = 0, (h) w = 0.6, (c) w = 1.2, and (d) 
w — 2.4. The energy gap is also color-coded such that black is zero gap and white is maximum on a scale from zero to 0.0425 
in Coulomb energy units, e^/^nl). 



duces the Coulomb energy by softening it. For w = 2.4, 
the maximum width considered, the energy gap retains 
the basic qualitative structure as for the other widths, 
however, its overall value is very small. Note as well that 
the maximum energy gap ridge is approximately located 
at the point in the d-t parameter space where the over- 
lap between the exact state crosses over from being larger 
with the Pfafhan ansatz and the 331 ansatz (more on this 
below), i.e., the gap is maximum close to the transition 
line between 331 and Pfaffian. 

In Fig. [5] we project the energy gap onto the two- 
dimensional t-d plane with the color coding indicating 
the numerical FQHE gap strength on the same scale as 
in Fig. 21 Note that the Coulomb interaction is undefined 
for d < w so no results are obtained in that case. On this 
same plot we draw with a black dashed line the line that 
separates the 331 phase from the Pfaffian phase. This 
phase boundary is determined by noting which state has 
a higher overlap with the exact state in which re- 
gion of parameter space, i.e., if (^331 |^o) is larger than 
(^Pf 1^0) then we say the exact system is in the 331 phase 
and vice- versa. Note that the dashed line is only an op- 
erational phase boundary within our calculation since all 
we know is that the 331 (Pf) has higher (lower) overlap 



above (below) this line. We cannot rule out the possi- 
bility that the in thermodynamic limit some other state 
(including possibly a compressible state) besides the ^331 
or ^pf dominates. 

We show the approximate quantum phase diagram 
(QPD) for all four values of the layer width parameter (a) 
w = 0, {h) w = 0.6, (c) w 1.2, and (d) w = 2.4. The 
zero- (Fig. [S^a)) and the intermediate- width (Fig. Efb)) 
results are of physical relevance whereas the (unrealisti- 
cally) large width results (Fig.[5fc) and (d)) are provided 
here only for completeness (since this is the regime where 
the Pf state dominates over the 331 state in the quan- 
tum phase diagram) . We note that we are using the sim- 
plistic Zhang-Das Sarma (ZDS) modelS for describing 
the well width effect, and crudely speaking w = 1 in the 
ZDS model corresponds roughly to wqw ^ 6 where wqw 
is the corresponding physical (i.e., effective single-layer) 
quantum well width. For a single WQW, where the ef- 
fective bilayer is created by the self-consistent potential 
of the electrons themselve^ii^i^i^, our w is typically 
much less than the total width W of the WQW-very 
roughly speaking w ~ and d ~ W/2. As empha- 

sized above, we treat t, and w{< d) as independent 
tuning parameters. 
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FIG. 5: (color online) Quantum phase diagram (QPD) and FQHE gap (color coded) versus layer separation d and tunneling 
strength t for widths (a) w = 0, (b) w = 0.6, (c) w = 1.2, and (d) w = 2.4. For the QPD, the 331 and Pf phases (as discussed 
in the text) are separated by a dashed black line and labeled appropriately. The FQHE gap is given as a contour plot with 
color coding given by the color-bar from dark to light, i.e., white being a largest value of 0.0425 and black being value of 0. 
The asterisks, triangles, circles, and squares correspond to the different experiments in Refs. fssl. Isol. l4ll. ITol and 111, respectively. 
Only experimental points showing FQHE are within the large solid circles in (a) with the lower smaller circles and upper larger 
circles indicating experiments in double-quantum-well structures and WQW structures, respectively. We note that the single 
triangle on the Pf side of the QPD does not manifest any experimental FQHE indicating that the theoretical gap may be 
overestimated for the Pf state. 



In Fig. \5\ we have put as discrete symbols all exist- 
ing u = 1/2 bilayer published FQHE experimental data 
(both for double quantum well systems and single wide- 
quantum- wells) in the literature, extracting the relevant 
parameter values (i.e., d and t) from the experimental 
works ^^^^-i^i^i^. Because of the ambiguity and uncer- 
tainty in the definition of w (i.e., how to precisely relate 
our theoretical w in the Zhang-Das Sarma model to the 
experimental layer width in real samples), we have put 
the data points on all four QPDs shown in Fig. [5] al- 
though the actual experimental width values correspond 
to only Figs. [5fa) and (b). See Sec. IVIIII for a detailed 
description of exactly how the experimental points were 
determined. 

Results shown in Fig. [5] bring out several important 
points of physics not clearly appreciated earlier in spite 
of a great deal of theoretical exact diagonalization work 
on = 1/2 bilayer FQHE: 

(i) It is obvious that large, or small, t and small, or large, 
d^ in general, lead to a decisive preference for the exis- 



tence of = 1/2 Pf, or 331, FQHE. The fact that large 
t values would preferentially lead to the Pf state over 
the 331 state is, of course, expected since the system be- 
comes an effective one-component system for large tun- 
neling strength. 

(ii) What is, however, not obvious, but apparent from 
the QPDs shown in Fig. [5l is that the FQHE gap (given 
in color coding in the figures) is maximum near the phase 
boundary between 331 and Pf. 

(iii) Another non-obvious result is the persistence of the 
331 state for very large (essentially arbitrarily large!) val- 
ues of the tunneling strength t as long as the layer sep- 
aration d is also large-thus having a large t by itself, as 
achieved in the Luhman et al. experiment^, is not enough 
to realize the single-layer v = 1/2 Pf FQHE, one must 
also have a relatively small value of layer separation d 
so that one is below the phase boundary (dashed line) 
in Fig. [5l The explanation for the Luhman experimental 
u = 1/2 FQHE being a 331 sate, as can be seen in Fig.[5l 
is indeed the fact that both t and d are large in these 
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FIG. 6: (color online) FQHE energy gap versus tunneling strength t for several values of layer separation d for (a) w = 0, (b) 
w = 0.6, (c) w = 1.2, and (d) w = 2.4. A dashed vertical line of the same color corresponds to the boundary between the 
Pfaffian phase (right of the line) and the 331 phase (left of the line). 



0.0.03 




FIG. 7: (color online) FQHE energy gap versus layer separation d for a few values of tunneling strength t for (a) it; = 0, (b) 
w = 0.6, (c) w = 1.2, and (d) w = 2.4. A dashed vertical corresponds to the phase boundary between the Pfaffian phase (left 
of the line) and the 331 phase (right of the line). 
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samples making 331 a good variational state. 

(iv) An important aspect of Fig. [5] is that the Pf FQHE 
gap tends to be very small-this is particularly true for 
larger values of where the Pf overlap is large. This 
implies, as emphasized by Storni et al?'^ ^ that the obser- 
vation of a = 1/2 Pf state is unlikely since the activa- 
tion gap would be extremely (perhaps even vanishingly) 
small. 

(v) For larger values of w (and large t), our calculated 
QPD is dominated by the Pf state-particularly for the 
unrealistically large width w = 2A (corresponding to 
"^QW ^ 14!) where all the experimental d and t val- 
ues fall in the Pf regime of the phase diagram. We em- 
phasize, however, that this Pf-dominated large-K; (and 
large-t) regime will be difficult to access experimentally 
since the FQHE gap would be likely extremely small as 
in Fig. [5l Our results however cannot decisively rule out 
the possibility of a = 1/2 Pfaffian state in the strong 
tunneling and small separation regime. 

We now discuss the published experimental results in 
light of our theoretical QPD. First, we note that most 
of the existing experimental points fall on the 331 side 
of the phase diagram which is consistent with our QPD 
in Fig. [5l In particular, only samples on the 331 side of 
the QPD with reasonably large FQH gaps, i.e., the data 
points close to the phase boundary, exhibit experimental 
FQHE. By contrast, the one data point (solid triangle in 
Figs.[5fa) and (b)) on the Pf side of the phase boundary 
does not manifest any observable FQHE in spite of its lo- 
cation being in a regime of reasonable FQHE excitation 
gap according to our phase diagram. This is consistent 
with the finding of Storni et al'^^ that the u = l/2 FQHE 
gap in a single-layer system is likely to be vanishingly 
small in the thermodynamic limit. It is, therefore, pos- 
sible that the Pf regime in our QPD has a much smaller 
excitation gap than what we obtain on the basis of our 
N = 8 particle diagonalization calculation. We refer to 
Storni et al^ for more details on the theoretical status 
of the single-layer LLL u = l/2 FQHE. 

For a more detailed view of the u = 1/2 bilayer FQHE, 
we show in Figs. [6l^a)-(d) and [7l^a)-(d), respectively, our 
calculated FQHE gap as a function of t (for a few fixed d 
values) and as a function of d (for a few fixed t values). 
Note that similar results were obtained by Nomura and 
Yoshioka in Ref. 18, however, they only considered the 
energy gap versus tunneling for two values of separation 
d and fixed w = 3.8 for di N = 6 electron system, which 
is aliased with a possible FQHE at = 2/3 and therefore 
suspect. In each figure, we also depict the line separating 
the 331 (smaller t/larger d) and the Pf (larger t/smaller 
d) regimes in the phase diagram. The qualitatively in- 
teresting point is, of course, the non-monotonicity in the 
FQHE gap as a function of t or d with a maximum close 
(but always on the 331 side) to the phase boundary. The 
non-monotonicity in the FQHE gap as a function of t (but 
not d) was earlier pointed out, but our finding that the 
peak lies always on the 331 side of the phase boundary is 
a new result. We emphasize that this result is strong ev- 



idence that the 331 phase is the dominant FQHE phase 
in u = 1/2 systems. We believe that the only chance 
of observing the u = 1/2 Pf FQHE is to look on the 
Pf side of the phase boundary at fairly large values of 
d and t. This is in sharp contrast to the SLL u = 5/2 
bilayer FQHE where we show in Section |V] that there 
are two sharp ridges far away from each other in the d-t 
space corresponding to the v = 5/2 Pf and 331 bilayer 
phases^. 

We conclude this section by commenting on the nature 
of the quantum phase transition (QPT) between the 331 
and the Pf phase in the d-t space. Our calculated QPD 
implies a continuous QPT from the strong-pairing 331 
to the weak-pairing Pf state with increasing t and/or de- 
creasing d as predicted by Read and Green^^. (Note this 
has also been very recently discussed by Dimov, Halperin, 
and Nayak^^.) In fact, our finite size diagonalization 
based QPD of Fig. [5] is topologically equivalent to the 
phase diagram predicted by Read and Green (Fig. 1 of 
Ref. 0). We emphasize, however, that we cannot distin- 
guish a quantum phase transition from a crossover be- 
cause of the limitations of finite size calculations. It is 
also possible that a different phase, e.g. a compressible 
composite fermion Fermi liquid phase, has lower energy 
and intervenes between the 331 and Pf phases so that 
the system goes from 331 to Pf (or vice versa) through 
two first-order transitions. What we have shown here 
is that if the u = 1/2 bilayer Pf phase exists at all, it 
would manifest most strongly in very wide samples and 
close to the phase boundary with the 331 phase with 
an extremely small FQHE excitation gap. We have also 
shown, through an explicit comparison with our t-d-w 
phase diagram (Fig.[5j) that all published = 1/2 FQHE 
dataiii^i^^i^Si^ are consistent with the existence of only 
the 331 phase in the LLL. Our work does not rule out 
the possibility of a weak Pf FQHE at z^ = 1/2 for large 
values of tunneling. 



V. SECOND LANDAU LEVEL 

We now consider the same set of questions we asked 
regarding the physics of FQHE bilayers in the = 1/2 
lowest LL (Sec. lIVp about the bilayer FQHE systems in 
the second Landau level (SLL). Our approach here is to 
theoretically consider a single-layer FQHE at a half-filled 
SLL system, i.e., = 2 + 1/2 = 5/2 where the first two 
lowest Landau levels of spin-up and spin-down are com- 
pletely occupied and inert. Thus, the half-filled SLL in- 
teracting electrons are projected into the LLL using the 
appropriate SLL Haldane pseudopotentials^^. The exact 
nature of this procedure has been given many times in 
many places in great detail and will not be reiterated here 
(see Ref. [30| for a good description). To consider a SLL 
bilayer system we then allow the single-layer to become a 
bilayer by reducing the tunneling and/or increasing the 
layer separation d. There are actually quite subtle points 
in defining this procedure in the SLL and they are dis- 
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FIG. 8: (Color online) Wavefunction overlap between the exact ground state and the Pf state ^pf (blue (dark gray)) and 
331 state ^331 (yellow (light gray)) as a function of distance d and tunneling amplitude t for the half- filled second LL with 
N = 8 electrons and width (a) w = 0, (b) w — 0.6, (c) w = 1.2, and (d) w — 2.4 (d > w necessarily). 




FIG. 9: (Color online) Wavefunction overlap between the exact ground state ^0 and ^pf at strong tunneling t = 0.2 (a) and 
^331 at zero tunneling t = (b) as a function of separation d and single-layer well width w (where d > w) for the SLL with 
N = 8 electrons. White corresponds to an overlap of unity while black corresponds to an overlap of zero. As in Fig. H the 
dashed red line is the condition w = d and for it; > d the bilayer system is undefined, i.e., the single- layer width cannot be 
larger than the layer separation. 

cussed below in Sec. IV Dl in detail. However, our pro- A. What is the physics?-second Landau level 

cedure is completely well defined theoretically, and is a 
direct analog of the u = 1/2 LL calculation in the SLL. 

The calculated wavefunction overlap between the ex- 
act ground state ^0 and the two appropriate candidate 
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FIG. 10: (Color online) {Ns — Na)/2 as a function of distance d and tunneling amplitude t for the half-filled second Landau 
level, N = 8 electrons, and (a) it; = 0, (h) w = 0.6, (c) w = 1.2, and (d) w — 2.4. 



variational wavefunctions (^pf and ^331) as a function 
of distance d and tunneling energy t is shown in Fig. [8] 
for single-layer widths of w = 0, w = 0.6, w = 1.2, and 
w = 2A. First we focus on the situation with zero width 
w = (Fig. Efa)) and concentrate on the overlap with 
^Pf. In the limit of zero tunneling and zero d the over- 
lap with ^pf is small-approximately 0.5 (unclear from 
the figure). However, only weak tunneling is required to 
produce a state with a sizable overlap of approximately 

0.96 and as tunneling increases this overlap remains 
large and approximately constant. For d 7^ 0, in the weak 
tunneling limit, the overlap with decreases drasti- 
cally. Adding moderate to strong tunneling we obtain 
a sizable overlap, decreasing gently in the large d limit. 
This shows that the strong-tunneling (one-component) 
regime is well described by the Pf state. This strong 
tunneling Pf regime appears to be robust. 

Next we consider the overlap between ^0 and ^331. In 
the zero tunneling limit, as a function of the overlap 
starts small, increases to a moderate maximum of ~ 0.80 
at d ^ 1 before achieving an essentially constant value of 
^ 0.56. For d > 4 the overlap remains relatively constant 
and slowly decreases as the tunneling is increased (in fact, 
there is a slight increase in the overlap to ~ 0.6 for a 
region of positive d > A and 0.1 ^ t ^ 0.15). Thus, the 
weak-tunneling (two-component) regime is well described 
by ^331. 

Similar to what we did in the LLL, we consider the 
effects of finite width, knowing that it may enhance the 
overlap^^'^^. In Fig. [9] (a) and (b) we show (^o|^33i) 



and (^ol^Pf), respectively, versus layer separation d and 
single layer width w for small tunneling t = and large 
tunneling t = 0.2. In the zero tunneling limit (Fig. [9fa)) 
we see that the overlap can be increased significantly 
(~ 0.94) by using d^ 1 and w ^ d. The overlap between 
the exact state and the 331 state remains high with in- 
creasing width as long as the layer separation d ~ 1-this 
compares with the LLL case where the overlap is maxi- 
mum for a larger separation near 1.8. However, the gen- 
eral feature that the maximum does not change apprecia- 
bly with increasing w is consistent with the LLL results. 
The main difference is the overlap with 331 is lower, and 
increasing w and d eventually pushes the overlap to zero. 

For the strong tunneling limit {t = 0.2), in the case 
of ^pf (Fig. El^b)), our results are consistent with the 
single- layer finite-thickness results (cf. Refs. 4Z andi25[). 
Compared to the LLL, however, the overlap increases to 
a maximum for a finite w and that behavior continues 
to hold when d is increased. Eventually, large single- 
layer width causes the overlap with the Pfaffian state to 
decrease for any d significantly larger than w. Thus, finite 
layer width enhances the overlap with both 331 and Pf 
states in their regimes of phase space. 

By examining Fig. [81(b)- (d) we see results qualitatively 
similar to the LLL results shown in Fig. [H The real dif- 
ference here is that the overlap between the exact state 
and the Pf state, in the region of phase space where it 
is better than the 331 state, is higher in the SLL than it 
is in the LLL. This is expected behavior when one con- 
siders the single-layer results^^^^^ where the Pf is known 
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FIG. 11: (Color online) FQHE energy gap as a function of layer separation d and tunneling amplitude t for the second LL 
system with N = 8 electrons. We consider single quantum well widths of (a) w = 0, (h) w = 0.6, (c) = 1.2, and (d) w = 2.4. 
The energy gap is also color-coded such that black is zero gap and white is maximum on a scale from zero to 0.035. 



to be an excellent candidate for the u = 5/2 single-layer 
FQHE. The other difference complimentary to the Pfaf- 
fian behavior is that the overlap between the exact state 
and the 331 state is lower in the SLL than it is in the 
LLL. In fact, for large layer width w = 2.4 (Fig. E^d)) 
the overlap with the 331 is quite low and it would be un- 
reasonable to assume that the exact state is adequately 
described by the Halperin 331 state. The Pfaffian overlap 
is also lower for w = 2A in the SLL than it is in the LLL 
(somewhat surprisingly). In fact, after investigating the 
FQHE energy gap for w = 2A (below in Sec. IV Cp it is 
clear that the bilayer SLL system most likely would not 
exhibit any FQHE for such large widths. Again, similar 
to the LLL results (Fig. [Hb)-(d)), we note that for in- 
creasing the region in phase space described by the Pf 
phase increases at the expense of the 331 state. 



B. Is the system one- or two-component?— second 
Landau level 

Fig. [To] shows the calculated value of {Ns — Na)/2 as 
a function of d and t with (a) w = 0^ {h) w = 0.6, (c) 
w = 1.2, and w = 2A and, as before, our overlap based 
conclusions are consistent with the value of (A^^ — A^^)/2. 



It is difficult to clearly tell the difference between the re- 
sults in the SLL compared to those in the LLL. The slight 
difference between the two is that in the large d limit 
slightly more tunneling is required to drive the system 
to the one-component regime than in the LLL case. Fur- 
ther, in the strong tunneling limit, the SLL system's one- 
component character is more robust to increasing layer 
separation d-but only just. Again, similar to the LLL, 
we find that when (A^^ — Na)/2 ~ 2.5 the overlap with 
the exact ground state switches from being either higher 
with the Pf state or the 331 state. 



C. Will the system display the FQHE?-second 
Landau level 

In Fig. [TlTa)-(d) we show the FQHE energy gap as a 
function of d and t for the SLL system for (a) = 0, (b) 
w = 0.6, (c) w = 1.2, and (d) w = 2.4. Similar to the 
LLL, it clear that for finite d and t, a 5/2 FQHE with 
a finite gap (being either 331 or Pf) exists in a realistic 
parameter regime. There is a clear qualitative difference 
between the lowest and second LLs, however. In the SLL, 
the largest FQHE energy gap is obtained in the zero d 
limit for finite tunneling. That is, in the region where the 
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FIG. 12: (color online) Quantum phase diagram (QPD) and FQHE gap (color coded) for the second LL for widths (a) w = 0, 
(b) w = 0.6, (c) w — 1.2, and (d) w — 2.4. For the QPD, the 331 and Pf phases are labeled appropriately and the gap is given 
as a contour plot with color coding given by the color-bar from dark to light, i.e., white being a largest value of 0.035 and black 
being value of 0. 



Moore-Read FfafEan ansatz is the better description of 
the exact ground state, the energy gap is largest. With 
an energy gap nearly as large is the "ridge" region identi- 
fied in our LLL results-again the ridge lies in the region 
of d-t phase space corresponding to the "quantum phase 
transition" between the Pf phase and 331 phase. Finite 
single-layer width w (Fig. [nT b)-(d)) decreases the overall 
energy gap and moves the ridge area to weaker tunneling. 
Note that for the largest width considered {w = 2.4) the 
energy gap is very small and the ridge has essentially be- 
come a "valley" and the overall energy gap is very small. 

We now discuss the approximate quantum phase di- 
agram for the bilayer system for the SLL as shown in 
Fig.[T2]for single-layer widths (a) il' = 0, (b) = 0.6, (c) 
w = 1.2, and (d) w = 2.4. This QPD is calculated the 
same way as it is for the LLL (cf. Fig. [5])- identifying the 
331 and Pf phases as the regions in the parameter space 
where the overlap with is larger (Fig. [8]). In general, 
the two-component 331 phase (the upper left region) has 
a weaker SLL FQHE than that of the one-component Pf 
(the lower right region) . Figure [121 which is an impor- 
tant prediction of our work, shows that in the SLL bilayer 
system, both the 331 and Pf states should be visible, and 
in a realistic finite thickness system (Fig. [T2fb) and (c)). 



the FQHE gap would become very small, perhaps even 
zero, in between the two phases. 

Again, we ask the natural question whether our nu- 
merically obtained transition from the large-d (small-t) 
region to the small-d (large-t) region is a quantum phase 
transition or a crossover. Of course a finite-size numer- 
ical study cannot definitively answer this question. Our 
results, however, are consistent with the findings^^^ of 
Read and Green, and we believe that there would be a 
QPT between 331 and Pf phases in the second LL since 
a crossover between topologically trivial (331) and non- 
trivial (Pf) phases is difficult to contemplate. We also 
note that for the w = case, the quantum phase tran- 
sition line seems to terminate at d = = t for the SLL 
and for d = and finite tunneling t in the LLL, i.e., some 
finite t is required to push the system into the Pf phase. 
Thus, our SLL results are more similar to the QPD of 
Read and Green^^ where a multi-critical point exists sep- 
arating the 331 Abelian phase from the non-Abelian Pf 
phase when t = 0. It is perhaps understandable that our 
SLL results would more closely resemble that of Read and 
Green since they were modeling a 1/2-filled bilayer FQHE 
state as Px + ipy superconductor of composite fermions 
which is thought to be the appropriate description for 
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FIG. 13: (color online) FQHE energy gap versus tunneling strength t for a few values of layer separation d for (a) it; = 0, (b) 
w = 0.6, (c) It; = 1.2, and (d) w = 2.4. A dashed vertical line of the same color corresponds to the boundary between the 
Pfaffian phase (right of the line) and the 331 phase (left of the line). 




FIG. 14: (color online) FQHE energy gap versus layer separation d for a few values of tunneling strength t for (a) it; = 0, (b) 
w = 0.6, (c) w = 1.2, and (d) w — 2.4. A dashed vertical line of the same color corresponds to the boundary between the 
Pfaffian phase (left of the line) and the 331 phase (right of the line). 
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the 1/2-filled FQHE in the SLL, not necessarily in the 
LLL. Further, in our SLL QPD, the Pf phase is quite 
strong-as indicated by the FQHE energy gap-along the 
d = and finite t line, whereas in the LLL, the gap is 
quite small along that line. Of course, the parameters in 
Read and Green's analysis are not directly related to our 
parameters-they have a chemical potential /i in an effec- 
tive pairing Hamiltonian while we have layer separation 
d. Furthermore, our model is SU(2) symmetric at the 
d = = t point while their model has explicitly broken 
SU(2). More work will be needed to understand com- 
pletely the similarities of our approach and that of Read 
and Green definitively, but the fact that the topology 
of our numerically calculated phase diagram of Fig. [121 
strongly resembles that of the phase diagram discussed 
by Read and Green is highly suggestive. 

In Figs. [T3r a)-(d) and [T4r a)-(d) we provide a detailed 
view of the SLL FQHE bilayer system with the FQHE 
energy gap given as a function of tunneling energy t for 
several values of the layer separation d and the energy 
gap as a function of layer separation for several values 
of tunnehng energy, respectively. The FQHE energy gap 
as a function of tunneling is similar qualitatively to the 
results in the LLL. Namely, the gap rises to a peak value 
and then falls off. Of course, for w = 2.4 (Fig. [TSTd)) 
the FQHE gap structure is different than all other cases, 
lowest or second LL, in that the peak has turned into a 
valley. 

In Fig. [14] we find that the behavior of the FQHE en- 
ergy gap as a function of layer separation d for constant 
tunneling energy t is qualitatively different from the be- 
havior in the LLL. The difference is that in the SLL the 
largest value of the energy gap is for d = and for increas- 
ing d the energy gap decreases slightly to a minimum (for 
t > 0.02) and then rises again to a peak. For increasing 
tunneling energy, the peak becomes more rounded. This 
behavior is qualitatively different from the lowest LL in 
that there are two peaks in the energy gap versus d at 
constant t instead of one-note that for w = and t < 0.02 
there is only a single peak not unlike the lowest LL re- 
sults. For zero layer width (k; = 0), the phase boundary 
between the one-component Pfaffian phase (to the left of 
the boundary line) and two- component 331 phase (to 
the right of the boundary) is to the left of the single 
peak in the energy gap for t < 0.02, moves to slightly 
right of the second energy gap peak for 0.02 < t < 0.1, 
and then finally moves back to being left of the second 
energy gap peak. For w = 0.6 in Fig. [T4r b) there are two 
energy gap peaks and the phase boundary is slightly to 
the right of the peak for t < 0.1 for w = 0.6. For t > 0.1 
the phase boundary moves to the left of the second peak. 
For w = 1.2 in Fig. [MT c) (except for t = 0) the phase 
boundary is always to the left of the second energy gap 
peak. Lastly, we show results for w = 2.4 in Fig. [T4lf d) 
for sake of completeness but it is clear that this case is 
most likely not a FQHE for any parameter values other 
than for very small d and large t. 



D. Bilayer FQHE in higher Landau levels 

There is a conceptual difficulty that we postponed 
when discussing bilayer FQHE in the SLL which exists 
when discussing any bilayer FQHE in higher LLs. In fact, 
we carefully set up the initial SLL bilayer FQHE problem 
as starting with a single-layer, or one-component, system 
in the SLL (presumably the 1/2 filled SLL at u = 5/2) 
and then systematically driving the system into a bilayer 
by the tuning of model parameters,. However, in our cal- 
culation, the electrons fractionally filling the SLL remain 
in the SLL throughout the procedure of tuning the sys- 
tem from one- to two-component. A bilayer FQHE at 
total filling factor u = 5/2 would consist of two layers 
at single layer filling factor z> = (1/2) (5/2) = 5/4, see 
Fig. [151 This is compared to the LLL where the half- 
filled bilayer FQHE consists of = 1/4 + 1/4 = 1/2, i.e., 
two layers at z> = 1/4, which when combined yield total 
iy = l/2 = 2v (Fig. [15Kb)). 

The case of spin-less electrons (we consider the spin- 
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FIG. 15: Schematic diagram depicting the FQHE system at 
total filling factor (a) z/ = 5/2 and (b) 1^—1/2 being broken 
up into a bilayer system with each separate layer being at 
(a) z> = 5/4 and (b) z> = 1/4. The colored shaded section 
represents a density of electrons filling up a quantum well. 
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FIG. 16: Schematic diagram depicting a spin-less FQHE sys- 
tem at total filling factor = 3/2 becoming a bilayer spin- less 
system with each layer at z> = 3/4. The fraction of the density 
fractionally filling a LL is colored blue (dark grey) while an 
inert LL (in this case the lowest) is colored red (light grey). 



full situation below) is depicted in Fig. [161 In the one- 
component limit, electrons in the half- filled SLL have to- 
tal filling factor v = 1 + 1/2 = 3/2-the LLL is completely 
filled and considered inert. In the two-component limit, 
when each of the two quantum well layers have been 
taken far apart from one another, the filling factor in 
each layer (assuming the electrons to have equal densi- 
ties in each layer) will be z> = 3/4. Thus, the electrons 
in individual layers are in the lowest LL, not the SLL. 
The key conceptual difficulty is that we started with a 
one-component system at total = 1/2 in the second LL 
and ended up splitting it into a bilayer system of each 
layer at z> = 3/4 in the lowest LL, i.e., the electron LL 
index changed! 

If one were to study bilayer FQHE in the SLL com- 
pletely rigorously (necessarily numerically) one would 
have to consider-even for spin-less electrons-a system 
consisting of electrons with a layer degree of freedom 
(pseudo-spin) and at least two LL degrees of freedom. 
Furthermore, one would not necessarily be able to treat 
the electrons in the LLL (the 1 in the 1 + 1/2 = 3/2) 
as inert. For the system we studied with = 8 elec- 
trons half- filling the second LL, one would need to fully 
consider N = 20 interacting electrons with a layer de- 
gree of freedom and consider at least two LLs. Hence, 
the Hilbert space dimension would be extremely large- 
on the order of 10^^ states. A Hilbert space dimension of 
10^^ is out of reach for any numerical procedure for any 
computer. 

The experimental reality makes the situation more 
comphcated due to the inclusion of spin. Half-filled 
FQHE in single-layer, presumably, one-component sys- 
tems has only been observed in the SLL, i.e., u = 5/2, 
which is modeled successfully by treating the completely 
filled spin- up and spin-down LLLs as inert. Once the 



system is made into a bilayer we are left with, in the ex- 
treme layer separation limit, two systems at 5/4 = 1 + 1/4 
filling. When the layers are very far apart (large d and 
necessarily weak t) it would be expected that the struc- 
ture of the 5/4 filled systems would be identical (since 
they would not interact) and the spin-up LLL would be 
filled and inert and the remaining electrons would fill the 
spin-down LLL up to 1/4 filling. Such a situation is de- 
picted in Fig. [TTTb). If, on the other hand, the system 
is barely becoming two-component, perhaps the simplest 
possibility is that the two z> = 5/4 layers will consist of a 
completely filled spin-up LLL in one layer and spin-down 
LL in the other layer (see Fig.fTTTa)). However, when the 
two layers are close enough that they begin to interact 
with one another-even without any tunneling-the situa- 
tion quickly becomes complicated (as discussed below). 
It is not clear if any of the electrons in a bilayer FQHE at 
total u = 5/2 can be considered inert and one might need 
to consider an interacting electron system where the elec- 
trons carry spin, pseudo-spin, and two, or possibly three, 
LLs indices, making such a problem intractable. 

We note that considering each layer to have z> = 2 + 1/4 
filling so that, in analogy with the LLL bilayer situa- 
tion, we have 1/4 filling in each SLL layer, does not re- 
solve the ambiguity since this produces a total filling of 

= (2 + 1/4) + (2 + 1/4) = 9/2, taking us to the half- 
filled third orbital LL! Similarly, considering each layer 
to have u = 1 + 1/4 filling leads to a total bilayer filling of 
v = 5/2, but takes us from the LLL to the SLL in going 
from a one layer to a bilayer system. The conundrum 
here is not theoretical, but is about how one connects 
the theoretical results with the experimental bilayer 5/2 
system. We emphasize that no such ambiguity arises for 
the = 1/2 = 1/4 + 1/4 bilayer system where all of the 
physics occurs in the LLL, both for individual layers and 
for the total bilayer. Also, there is no ambiguity in the 
strong tunneling regime where the bilayer effectively acts 
as a single layer system with z/ = 5/2 as in a single 2D 
system. 

For our fully spin-polarized model, however, there is an 
easy conceptual (and operational) way out of this ambi- 
guity for the bilayer u = 5/2 system. Since we assume the 
whole system (in both layers and in both orbital Landau 
levels) to have the same spin, the system is effectively 
spin-less (which is equivalent to assuming the Zeeman 
energy to be much larger than the cyclotron energy). In 
this spin-less (or spin-polarized) situation, each orbital 
LL by definition comes with just one spin index. There- 
fore, the V — 5/2 balanced bilayer system is equivalent 
to a z> = (1/2) (5/2) = 5/4 = 1 + 1/4 filling in each layer, 
where the completely filled inert level in each layer is sim- 
ply the LLL, and the 1/4-filled level in each layer is in 
the SLL. Since the LLL is inert, the incompressible FQH 
states form entirely in the SLL, and we can construct 
a Halperin 331 SLL state with v — \/2 from the two 
1 / 4-filled SLL states in each individual layer in a manner 
similar to that for the = 1/2 LLL bilayer state. We 
emphasize that this simplicity is lost if we include both 
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We mention that earlier theoretical work (alluded to 
above) by Zheng et al.^^ ^ Das Sarma et a/.,^^'^^, Demler 
et al.^^ ^ and Brey et al.^^ did take into account the dy- 
namical interplay between layer and spin indices with the 
conclusion that there should be a novel quantum canted 
antiferromagnetic phase in bilayers for situations where 
each individual layer state has a gap in the spectrum 
(e.g., z/ = 1 + 1 or 1/3 + 1/3 or 1 + 1/3 + 1 + 1/3). For the 
case we are considering in the current work, i.e., v = b/2 
and V = u/2 = 5/4, there is not necessarily a gap in the 
individual layer spectrum, and therefore a canted phase 
is not expected. 



(a) I d I 



VI. CONCLUSIONS 
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FIG. 17: Schematic diagram depicting two possible scenarios 
for the spin-full FQHE system at total filling factor = 5/2 
being broken up into a bilayer system with each separate layer 
being at i> = 5/4. In (a) we consider the small d and strong t 
limit where the lowest spin-up and spin-down Landau levels 
are completely occupied leaving a 1/4 filled spin-up Landau 
level in each layer. In (b) we depict the large d and weak 
tunneling t limit with the lowest spin-up LLL being filled in 
each layer leaving a 1/4- filled lowest Landau level of spin- 
down electrons in each layer. The notation is such that T^((j,n) 
corresponds to the filling factor i/ for electrons with spin-cr 
electrons with orbital Landau level index n, i.e., 1(^,0) corre- 
sponds to filling factor 1 of spin-t electrons with LLL index 
0. 



spin and layer indices on an equal footing since includ- 
ing the two spin indices (up and down) and two layer 
indices (right and left) lead to an immediate problem on 
how to assign individual z> values which, when combined 
into the total u = 2z>, lead to a bilayer SLL u = 5/2 
FQHE with both the total bilayer (i.e., u) and the in- 
dividual layer (i.e., z> = u/2) filling factor being in the 
SLL. (As emphasized already, the situation with two lay- 
ers, two spins, and two orbital LLs is ambiguous unless 
the numerical diagonalization takes exactly into account 
the full dynamics of all three two-level quantum indices 
which is impossible to do for any computer.) 



Our results concerning the u = 1/2 LLL FQHE bilayer 
system point to strong evidence that the dominant FQH 
phase is the two-component Abelian Halperin 331 phase. 
We believe that the only chance of observing the = 1 /2 
non-Abelian Moore-Read Pfaffian FQHE is to look on 
the Pf side of the phase boundary (Fig. [5]) at fairly large 
values of d and t. This contrasts the SLL bilayer situation 
where there are two sharp ridges far away from each other 
in the d-t space corresponding to the Pf and 331 phases. 
We note that for unrealistically large single layer width 
w the Pf phase dominates the 331 phase but the FQHE 
gap is extremely small. If the Pf phase exists at all, it 
would manifest most strongly in wide samples and close 
to the phase boundary with the 331 phase. 

In addition, we predict the existence of both the two- 
component 331 Abelian (at intermediate to large d and 
small t) and the one-component Pfaffian non-Abelian (at 
small d and intermediate to large t) u = 5/2 SLL FQHE 
phase in bilayer structures. The observation of these two 
topologically distinct phases, one (Pf) stabilized by large 
inter-layer tunneling and the other (331) stabilized by 
large inter-layer separation would be a spectacular verifi- 
cation of the theoretical expectation^^ that bilayer struc- 
tures allow quantum phase transitions between topologi- 
cally trivial and non-trivial paired even-denominator in- 
compressible FQHE states. The direct experimental ob- 
servation of our predicted "two distinct branches" of two 
strong SLL FQHE regimes in bilayer structures, as shown 
in Fig. [T2fa)-(c), with the FQHE gap being largest along 
the two ridges in the phase diagram as d and t are var- 
ied, will be compelling evidence for the existence of the 
u = 5/2 non-Abelian Pf state. 

Note that for small d and t, Pf and 331 phases are con- 
tinuously connected, indicating the possible existence of a 
quantum phase transition. We cannot, of course, dismiss 
a crossover due to the limitations of finite size calcula- 
tions, but it is difficult to contemplate how an Abelian 
and non-Abelian phase could continuously go into each 
other without a quantum phase transition. 

Besides our exact diagonalization results for the bilayer 
u = 1/2 and 5/2 FQHE, we have raised an important con- 
ceptual issue involving the existence of Halperin type bi- 
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layer FQHE states in higher (i.e., beyond the lowest) Lan- 
dau levels. In particular, the fact that two well-separated 
distinct layers have individual filling factors which 
necessarily lie in a lower orbital LL, make it tricky to 
define a composite bilayer Halperin wavefunction for v 
which is in a higher LL. (This problem obviously does 
not arise in the LLL.) It is conceivable that any theoreti- 
cal study of bilayer FQHE states in higher LLs must nec- 
essarily include the full dynamics involving layer, orbital 
LL, and spin degrees of freedom. This is an almost im- 
possible theoretical challenge, and we hope that our work 
will motivate experimental activity in bilayer structures 
at total filling factor 5/2 in order to explore the possi- 
ble difference between the physics of bilayer = 1 /2 and 
= 5/2 systems. 
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VIII. APPENDIX: IDENTIFYING THE 
EXPERIMENTAL POINTS IN THE zy = 1/2 
BILAYER FQHE QUANTUM PHASE DIAGRAM 

We briefly describe our procedure for identifying the 
experimental data points in our calculated = 1/2 bi- 
layer FQHE quantum phase diagram (Fig. [5]). In other 
words, we explain how we extracted the parameter val- 
ues t (tunneling energy), d (layer separation), and w 
(layer width) for each experimental sample point shown 
in Fig. [5l 

The tunneling energy t is connected directly to the 
symmetric- antisymmetric gap A^^^ = t in our model 
and therefore we simply use the quoted value of t (or 
^SAs) from the relevant experimental paper. The layer 
separation d appears in our int erlayer Coulomb interac- 
tion as l^nter = 6^/ + d^). For a true bilayer sys- 
tem the definition of d is obvious (the well-center to well- 
center distance between the two quantum wells), and we 
use that value for d from the experimental publications 
(see Fig. [TSTa)). For a wide single- well system, where 
the effective bilayer arises as a direct result of the self- 
consistent density profile in the system, we choose d as 
the distance between the two peaks in the calculated car- 
rier density profile for the corresponding wide well sam- 
ple (see Fig.fTSTb)). These self-consistent density profiles 
are given in each experimental paper for each sample we 
show in Fig.O This procedure is, of course, imperfect-for 
example, the density profiles could be different from the 
calculated zero-field local density approximation (LDA) 
results since the systems is in the strong field FQHE 
regime. We do not, however, believe that this is a serious 
issue since the 2D dynamics of the electrons leading to 
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FIG. 18: Schematic diagram depicting the two bilayer struc- 
tures considered in this work, (a) shows a double quantum 
well structure along with a generic density profile in the z- 
direction n[^z) (blue dashed line) for experiments like Ref. IH 
while (b) shows a single wide-quantum-well structure and 
generic density profile (blue dashed line) typical of experi- 
ments like Ref. [39l. 



Experiment 


d/l t/{ 


eV(^O) 


w/l 


Luhman et aL, (squares), Ref. ^ 


5.3 


0.13 


0.4 




6.9 


0.07 


0.3 


Eisenstein et aL, (asterisks), Ref. 38 2.4 


0.01 


0.3 




2.7 


0.01 


0.3 




3.6 


0.01 


0.3 


Suen et aL, (triangles), Ref. 39.41 


5.3 


0.08 


0.4 




4.8 


0.09 


0.4 




3.9 


0.14 


0.4 


Shabani et a/., (circles), Ref. 19 


5.9 


0.11 


0.25 




6.0 


0.13 


0.25 



TABLE I: Values of layer separation quantum well width 
w, and tunneling energy t in units of magnetic length / and 
Coulomb ener gy respectively, taken from Refs. [l|, 

m, S [U, andlii as described in Section Ivml 



FQHE is separable from the calculated density profile in 
the z-direction transverse to the 2D plane. In any case, 
the value of as given by the calculated LDA density 
profile, is the best one can do. 

The most problematic parameter for our model is the 
finite width parameter w which, in the Zhang-Das Sarma 
model (Ref. 57), does not correspond at all to the phys- 
ical layer width. This is why we have shown the exper- 
imental points on all four "width" values in Fig. [5l It 
is known from earlier work^« i^^i^^~ — that the Zhang-Das 
Sarma parameter w has the following, very approximate. 
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correspondence with the quantum well width parameter 
'^QW- — ^ ' Given that a balanced, symmetric, 
wide-quantum-well structure with a total width of W 
has wqw = W/2 for each individual layer, we conclude 
'^ZDS ^ '^Qw/12. This is, however, a very qualitative 
and crude estimate. Because of all the approximations 
involved in our depiction of the experimental data, points 
in Fig. [5] should be taken as a qualitative comparison, 
rather than a quantitative one. 

In Table H] we provide the values for t, and w we get 
from individual experimental papers whose data points 
show up in Fig. [5] 



Lastly, we schematically show in Fig. [18] how d and 
w in our model correspond to those in the experimental 
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